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ABSTRACT 

Summary: BioNetGen is an open-source software package for 
rule-based modeling of complex biochemical systems. Version 2.2 
of the software introduces numerous new features for both model 
specification and simulation. Here, we report on these additions, 
discussing how they facilitate the construction, simulation, and 
analysis of larger and more complex models than previously possible. 
Availability: Stable BioNetGen releases (Linux, Mac OS/X, and 
Windows), with documentation, are available at http://bionetgen.org 
Source code is available at http://github.com/RuleWorld/bionetgen 
Contact: bionetgen.help@gmail.com 

Supplementary information: Supplementary materials are available 
at Bioinformatics online. 


1 INTRODUCTION 

Rule-based modeling is an approach for addr essing combinatorial 
comp lexity in models of biochemical systems l lChvlek et all 1201 4 
I2OI5I) . Instead of manually enumerating all possible species 
and reactions that can exist within a system, a rule-based 
model defines only the reactive motifs within macromolecular 
complexes and the interactions and modifications that involve 
those motifs. BioNetGen is an open-source soft ware package 
for c o nstructing and simu lating rule-based models l lBlinov et ali 
|2004 iFaeder et ali l2009h . It has been used to model a variety 
of biological pr ocesses, including cell sig naling, gene regulation, 
and metabolism dChvlek et al[ 12014 1201 4 and references therein; 
see also bionetgen.org/index.php/ModeLExamplesl. Models are 
written in a human-readable, text-based modeling language known 
as BNGL (BioNetGen language). Numerous user-specified actions 
can be added to a BNGL model file, including generating 
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Fig. 1. Basic overview of the RuleBender/BioNetGen/NFsim software 
stack. BNGL models may be constructed and edited in RuleBender as well as 
translated from an SBML model using Atomizer (either with the BioNetGen 
readFile action or the web application atratomizer.appspot.com 1 . Models 
can then he simulated in three ways: (1) with a built-in population-based 
network simulator (ODE, SSA, PEA) after explicit network generation; (2) 
as a hybrid particle/population (HPP) model in NE sim after pa r tial ne twork 
expansion (“lumping” species must be provided iHogg et ali l2014h l: (3) 
with NEsim directly. In all cases, output trajectories are passed back to 
RuleBender for display. BioNetGen also has a number of built-in methods 
for exporting to third-party formats, such as SBML and Matlab. (ODE: 
ordinary differential equations; SSA: stochastic simulation algorithm; PLA: 
partitioned-leaping algorithm; GML: Graph Modeling Language.) 


a reaction network and performing deterministic or stochastic 
simulations. Mod els can also be exp orted to different formats, 
such as SBML dHucka et a/.L 1200.31 ) and MATLAB language 
(The MathWorks Inc., Natick, MA). Furthe rmore, BioNetGen 
interfaces with NEsim dSneddon et ali |20I ih . a “network-free” 
simulator that avoids enumeration of species and reactions, 
which may be intractable for large models. RuleBender is a 
visual interface for BioNetGen, with features that include syntax 
checking/hi ghlighting and visualizations for model de bugging and 
comparison dXu et al. J201 iLlWenskovitch et fl/.L|20l4 . BioNetGen 
is also used as a network generator and sim ulator in a number o f 
third-party tools including the V irtual Cell {Marwnet^, 2008), 
BioUML dKolpakov et ali 1 20061) . SRSim GrUnertandDittri^ 
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I 2 OI ih . Parts & Pools ( MMrhisi^ 2014h . pySB llLopez et fl/.Ll2013h . 
and BioNetFit llThomas et al. [ I 2 OI 6 I) . 

Taken together, the BioNetGen/NFsim/RuleBender suite of 
tools (Fig.lTJ provides powerful capabilities for rule-based modeling 
and simulation. Until recently, however, there were several 
significant shortcomings: (i) BNGL models were limited to 
mass-action kinetics; (ii) network-based stochastic simulations 
were limited to computationally expensive “exact” methods; 
(iii) network-free simulations were limited to fewer than a 
few million particles due to high memory load; (iv) there 
was no facility for importing models from SBML format. 
Each of these shortcomings has been addressed in the recent 
BioNetGen 2.2.x series of releases, as described below. Additional 
information can be found in extensive online documentation at 
http://bionetgen.org/index.php/Documentation (see Supplementary 
Information for specific links). 


2 ADDITIONS 

2.1 Functional rate laws 

BioNetGen 2.2 introduces the ability to define arbitrary mathematical 
functions that can be used as rate laws. Functional rate laws 
are integrated into the BNGL grammar and are supported by 
both NFsim and the network-based simulators (deterministic and 
stochastic) included in the BioNetGen package. Functions are 
defined using “observables,” wh ich compute the co ncentrations of 
species with specified properties l lFaeder et fl/.Ll200^ . and can either 
be evaluated globally, i.e., over the entire system, or locally, over a 
specified molecule or molecular complex. Local functions greatly 
expand the expressiveness of BNGL because they enable a single 
rule to specify many reactions wit h rates that depend on the specific 
properties of the reacting species iSneddon et fl/.Ll2oTll) . 

2.2 Accelerated-stochastic simulation 

T leaping is an approach for accele r ating stochastic simulations 
of biochemical networks dGillesniel l2007h . Numerous variants 
of r leaping have been proposed in the literature, including a 
multiscale variant known as the partitioned-leaping algorithm (PLA) 
^Harris and Clanc\ll200fih . In their simplest realizations, r-leaping 
methods are analagous to the explicit forward Euler method for 
solving ordinary differential equations (ODEs). Therefore, as with 
ODE integrators, higher-order and implicit versions of r leaping 
algorithms are possible. BioNetGen 2.2 includes an explicit Runge- 
Kutta implementation of the PLA (RK-PLA) that can be used on any 
model for which a reaction network can be generated. Additional 
information, including a performance analysis, is provided in the 
Supplementary Information. 

2.3 Hybrid particle/population simulation 

Network-free simulation methods that do not enumerate species 
and reactions are often required to simulate complex models 
iChvlek et a/.L I20l3 . l2015l) . However, they are limited by the fact 
that memory usag e increases linearly with the number of particles. 
iHogg et ^ ( I 2 OI 4 I) introduced a hybrid simulation method that treats 
a user-defined set of species as population variables rather than as 
particles. This hybrid particle/population (HPP) approach avoids 
the memory costs associated with having large pools of identical 


particles and was shown to significantly reduce computational 
memory expense with no effect on s imulation accuracy and little 
effect on run time jHogg et a/.Ll2014l) . HPP is implemented within 
BioNetGen 2.2 and can be run using NFsim version 1.11 or later. 


2.4 SBML-to-BNGL translation 

SBML i s a widely-used model exchange format in systems 
biology jHucka et a/.L l2003l) . Models encoded in SBML are flat, 
i.e., their species do not have internal structure, which limits 
their utility for rule-based modeling. BioNetGen 2.2 includes an 
SBML-to-BNGL translator, called Atomizer (also available as 
a web tool at ratomizer.appspot.com), that can extract i mplic it 
molecular structure from flat species l lTapia and Faedeii l2013l) . 
A full re port on Atomiz er and its application to the BioModels 
database i LLeLfij, l2010h is currently in preparation. However, 
iTapia andFaete ( 20131) reported that an early version of the tool 


could recover implicit structure for about 60% of species in models 
within the database that contain >20 species. Thus, Atomizer 
makes available a large set of pre-existi ng models in a rule-based 
format, facilit ating their visualiza tion llWenskovitch et al\. I20l4l 
and extension dChvlek et fl/.Ll201^ . 


2.5 Additional features 

BioNetGen 2.2 also introduces a number of additional features 
to those described above, including a null symbol (‘0’) that 
can act as a source or a sink in rules and new actions for 
performing parameter scans, generating MATLAB Executable 
(MEX) files, exporti ng to formats readable by third-party 
simulators (e.g., MCell l lKerr et aZ , 2008|)). an d outp utting graphical 
visualizations at different scales Isekar et fl/ll20lk) . We have also 
developed parameter estimation tool called ptempest (available 
at http://github.com/Ru leWorld / ptemp e stl that combi n es Ba yesian 
Monte Carlo meth odsi 


Monte Garlo meth odsaKlmkel. IzUI^; IEy( 
parallel tempering llEarl and PeemL 200511 


leWorld / ptemp e stI that combi n es Ba yesian 
Klinkel. boiA lEydgahi et all 1201 3l) with 


to accelerate sampling. 
Many new simulation options have also been added, e.g., to allow 
simulation arguments to be loaded from file and to terminate 
a simulation if a user-defined logical condition is met. A 
comprehensive listing of all new actions and arguments added in 
BioNetGen 2.2 is provided in the Supplementary Information (see 
also http://bionetgen.org/index.php/BNG_Actions_Args I. 


3 CONCLUSION 

BioNetGen is an active, open-source project that highly encourages 
contributions from interested members of the modeling community. 
Contact information, as well as numerous links to online 
documentation, is provided in the Supplementary Information. 
Ongoing d evelopment eff o rts in clude interfacing with spatial 
simulators |Sullivan et g/.L l2015h . free-energy-based modeling 
llHoggll20L3h . and improved support for community standards (e.g., 
SBML-multi, SED-ML). Recent independent efforts that leverage 
the BioNetGen framework include a standard for annotation of rule- 
based model s ( Mjsirli et g/.L 1201 6h specified in BNGL or Kappa 
dPanos et g/.L |2007|) . anoth er widely-used rule-b ased modeling 
language, and BioNetFit dThomas et al. I Iroleh , a parameter 
estimation tool for models simulated using either BioNetGen or 
NFsim. 
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